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00 ■ Abstract 



Kaluza-Klein (KK) black hole solutions in the Einstein-Gauss-Bonnet (EGB) grav- 
ity in D dimensions obtained in the current series of the works by Maeda, Dadhich 



$_i ' and Molina are examined. Interpreting their solutions, the authors claim that the 

mass/matter is created by the extra dimensions. To support this claim, one needs to 
■ show that such objects have classically defined masses. We calculate the mass and 

' mass flux for 3D KK black holes in 6D EGB gravity whose properties are sufficiently 



physically interesting. Superpotentials for arbitrary types of perturbations on arbi- 
trary curved backgrounds, recently obtained by the author, are used, and acceptable 
mass and mass flux are obtained. A possibility of considering the KK created matter 
as dark matter in the Universe is discussed. 

1 Introduction 

We study new exact solutions in the Einstein-Gauss-Bonnet (EGB) gravity in D dimensions, 
which are <i-dimensional Kaluza-Klein (KK) black holes (BHs) with (D — <i)-dimensional 
submanifold, presented recently in [1] - [4] by Maeda, Dadhich and Molina. The authors 
treat them as a classical example of creating matter by curvature. The idea of such a kind is 
not new. Thus, to make inflation possible, a pioneer proposal was advanced by Starobinsky 
[5] that a high-energy density state was achieved by curved space corrections. Many other 
problems of modern cosmology may be solved in the framework of multidimensional gravity 
using high-order curvature invariants of KK type spacetimes, see, e.g., [6] and references 
there in. 

To support the claim on creating 'matter without matter', it is necessary to calculate 
the mass and the mass flux by classical methods. It is the main goal of the present paper. 
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Here, we concentrate on 3D BHs in 6D EGB gravity [4]. These toy objects are rich enough 
in physical properties, e.g., they can have a radiative regime. For calculations we use the 
conservation laws developed by us in [7] - [9] , where in the framework of EGB gravity, super- 
potentials (antisymmetric tensor densities) for arbitrary types of perturbations on arbitrary 
curved backgrounds have been constructed. Three important types of superpotentials [9] are 
used, those based on (i) Ncether's canonical theorem, (ii) Belinfante's symmetrization rule 
and (iii) a field-theoretical derivation. 

The paper is organized as follows. In section 2, we outline the solutions obtained in [1] - 
[4] and describe necessary properties of the 3D objects in 6D EGB gravity. In particular, in 
a natural way, we define a spacetime where a BH is placed. It can be considered as a possible 
background against which perturbations are studied. In section 3, in the preliminaries, the 
main notions and properties of the applied formalism are presented. Then we study the 
objects themselves: (a) as vacuum 6D solutions; (b) as 3D KK solutions with a 'matter' 
created by extra dimensions. Calculating the mass and the mass flux we support the second 
viewpoint. In section 4, we discuss (a) an ambiguity in the canonical approach related to a 
divergence in the Lagrangian; (b) a possibility of applying the KK BH solutions in cosmology. 
The Appendix presents explicit general expressions for all three types of superpotentials in 
EGB gravity. 



2 Kaluza-Klein 3D black holes 



We consider the action of the EGB gravity in the form: 



S = / d D xt EGB = - J- / d D x, 
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R -2A + a (R^-AR^ + R 2 ) 



(2.1) 



where a > 0. Here and below, curvature tensor R^ vprj) Ricci tensor R^ and scalar curvature 
R are related to the dynamic metric g pu ; a 'hat' means densities of the +1, e.g., g^ u = 
V~99^j = d& means ordinary derivatives; the subscripts and 'gb' are related to 
the Einstein and the Gauss-Bonnet parts in (2.1). 

The main assumption in [1] - [4] is that the spacetime is locally homeomorphic to A4 d x 
K, D ~ d with the metric g^ v = diag((?AB, r olab), A,B = 0, •••,d— 1; a, b = d, • • • , D — 1. 
Thus, qab is an arbitrary Lorentzian metric on Ai d , 7 a & is the unit metric on the (D — d)- 
dimensional space of constant curvature K D ~ d with k — 0, ±1. Factor r is a small scale of 
extra dimensions compactified by appropriate identifications. The gravitational equations 
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corresponding to the EGB gravity action (2.1) have the form: 



Q\ = G" v + olH\ + <^A = , 



(2.2) 



where the Einstein tensor G^y and 5^ u correspond to the Einstein part and if ^ corresponds 
to the GB part in (2.1). After all assumptions their decomposition is as follows: 
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5 A b = 0; (2.3) 



(D-d-l)(D-d-2)x 







(2.4) 



where the subscript '(d)' means that a quantity is constructed with the use of Qab only. As a 
result, one can see that (2.3) is a tensorial equation on M d , whereas (2.4) is a constraint for 
it. However to obtain more interesting solutions one has to consider a special case that the 
quantity Q A b disappears identically. This is possible for d < 4 only because then (d)Hnv = 0. 
Next, constants are chosen so as to suppress the coefficients in (2.3), which is possible if 
D > d + 2, k = — 1 and A < 0. Taking into account all the above, there remains a single 
governing equation, the scalar equation (2.4) on M d . 

Here, we consider the solutions for D = 6 and d = 3 presented in [4]. A suitable set 
of constraints for the constants is Tq = 12a = — 3/A . Then, the left hand side of (2.3) 
disappears identically. Keeping in mind that (z)Lgb = 0, one simplifies (2.4) to obtain 



(d)R = 2A , 

to which the static solution <?ab(t") has been found: 

ds 2 = -fdt 2 + f~ l dr 2 + r 2 d<j) , / ee r 2 /l 2 + q/r-fi. 



(2.5) 



(2.6) 



Here, \x and q are integration constants, and I 2 = — 3/A . The Einstein tensor components 
for the solution (2.6) are 



G° = G\ = l/l 2 - q/2r\ (% = l/l* + q y 



(2.7) 



As a space of a constant curvature, (D — d = 3)-sector is completely presented by its scalar 
curvature: 

{D _ d) R = 6k/r 2 = 2A = -l/2a . (2.8) 
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For comparison we consider the BTZ BH [10]. Its metric is presented in the form 

ds 2 = -fdt 2 + f- x dr 2 + r 2 d(j), f = -r 2 A -/x, (2.9) 

which is a solution to the 3D pure Einstein equations. The horizon radius r + of the BH 
is defined as r\ = —fx/A , thus r + (and consequently a BH itself) disappears for vanishing 
pi. Therefore the integration constant /x can be called the mass parameter. For /x — > 0, the 
so-called real vacuum related to the BH (in another word, a spacetime where a BH is placed) 
is defined by (2.9) with / = — r 2 A . However, such a spacetime is not maximally symmetric, 
unlike AdS one. The latter with / = — r 2 A + 1 is approached when fx = — 1. A difference 
between a real vacuum and a maximally symmetric vacuum is usual in BH solutions of 
modified metric theories (see, e.g., [11, 12]); the BTZ BH is the simplest illustration. 

The solution (2.6) is more complicated than (2.9), although one has clear analogies with 
the BTZ case. Considering BH solutions for simulating dark matter (see a discussion in 
section 4) we are more interested in the cases with a horizon. In (2.6), the equation for the 
event horizon is l 2 q + r + (r+ — I 2 fx) = 0. It is again natural to choose a mass parameter fx 
in such a way that the BH horizon disappears under vanishing fx. This gives fx — fx — q/r + 
and r\ = l 2 Ji (compare with the BTZ case), and consequently fx > 0. Then a real vacuum 
is defined by (2.6) with / = r 2 /l 2 + q/r — q/r + , it is again not maximally symmetric. The 
maximally symmetric AdS vacuum is defined by (2.6) with / = r 2 /l 2 + 1. For the latter, 
parameter q is considered entirely as a perturbation together with fx + 1. For fx < a horizon 
does not exist, this takes place, when fx > with q > 21 (/i/3) 3 ^ 2 or /x < with q > 0. 

The scalar equation (2.5) is also satisfied by the radiative Vaidya metric gAB(v,r): 

ds 2 = -fdv 2 + 2dvdr + r 2 d<p, f = r 2 /l 2 + q(v)/r - fx(v) (2.10) 

where /x(v) and q(v) now depend on the retarded/advanced time v. Keeping in mind a 
possibility to form KK black holes [1] - [4], advanced time is more interesting. Then (2.10) 
can be connected with the solution of the form (2.6) by the transformation dt = dv — 
dr/f(v,r). After that, for every constant vq, one can define its own horizon (if it exists) and 
a corresponding real vacuum analogously to the static case. The Einstein tensor components 
corresponding to (2.10) are 

G° = G\ = \jl 2 - q/2r 3 , Gl = {J^r-q)/2r 2 , G\ = l/l 2 + q/r 3 , (2.11) 

where dot means d/dv. The scalar curvature of (D — d = 3)-sector is expressed again by 
(2.8). 

Considering (2.6) and (2.10) as solutions to the Einstein 3D equations on M 3 (or, the 
same, EGB equations because in (2.2) one has (^H^v = 0), one concludes that they are not 
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vacuum equations with a redefined cosmo logical constant A = A /3 = —I /I 2 . Indeed, both 
(2.7) and (2.11) show that a 'matter' source Tab with zero trace T A a = should exist, and 
the Einstein equations corresponding to (2.5) could be rewritten as 

(3)Rab — Iqab(3)R + 5UbA = k 3 Tab ■ (2-12) 
A natural treating in [1] - [4] is that Tab is created by the compact extra dimensions. 

3 The mass and the mass flux for 3D black holes 
3.1 Preliminaries 

Our calculation is based on differential conservation laws for perturbations in a given back- 
ground spacetime in the form: 

l a (0=d^(0 (3.1) 

where £ a is a displacement vector, T a is a vector density (carrent) and T a/3 is an antisymmet- 
ric tensor density (superpotential). Thus, d a pT a ^ = and d a T a = 0. The current contains 
energy-momentum of both matter and metric perturbations, whereas the superpotential de- 
pends on metric perturbations only. Integrating d a T a = and using the Gauss theorem one 
obtains the integral conserved charges in a generalized form: 

V(0 = I d^xF® = I dSi$*(t) (3.2) 

JdT, 

where £ is a (D — 1) -dimensional hypersurface x° = const, <9£ is its (D — 2)-dimensional 
boundary, the zero indices correspond to time or lightlike coordinates, and small Latin indices 
correspond to space coordinates. Since we consider spherically symmetric systems, we need 
01-components of the superpotentials in (3.2) only. 

The formalism describes exact (not infinitesimal) perturbations in general. This is 
achieved if one one solution (dynamical) is considered as a perturbed system with respect 
to another (background) solution of the same theory. Thus conserved quantities are defined 
with respect to a fixed (thought as known) spacetime, e.g., a mass of a perturbed system 
on a given background. A background can be both vacuum and non-vacuum, and usually 
is to be chosen to correspond with problems under consideration. The task of the present 
paper is calculating a global mass of the KK BHs presented above. It is more important 
the mass defined with respect to a spacetime, in which BH is placed because then with 
vanishing BH, one obtains a zero mass. Therefore, first of all a real vacuum described in 
previous section is chosen as a natural background. Although such backgrounds are curved 
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and nonsymmetric, the technique used is powerful. Besides, as interesting and important 
backgrounds we consider the AdS space. For such kinds of backgrounds, perturbations are 
not infinitesimal in general. However, we need in appropriate asymptotic of superpotentials 
in (3.2) only. As one can see below, the fall-off integrands in (3.2) both at spatial and at 
null infinity turns out to be sufficiently strong to allow surface integrals to converge and to 
give reasonable results. 

In the previous section, the bar meant a quantity related to a spacetime where a BH 
is 'placed'; here and below, without contradictions the bar means a quantity related to a 
background spacetime as a structure of the formalism. As a natural choice, for the above 
described static and radiative solutions we use the background metric in the same forms 
(2.6) and (2.10), respectively, where / = f(r) can be arbitrary in general but should be 
static. For calculating the global mass M we use the timelike Killing vector 

r = (-l,0). (3.3) 

It has this unique form for the above two generalized types of background metrics: the zero 
component in (3.3) can be both timelike and lightlike; includes 5 or 2 space dimensions 
in a 6D or 3D derivation, respectively. The metrics of the real vacuum and AdS space just 
belong to the aforementioned two types of background metrics and consequently also have 
a timelike Killing vector of the unique form (3.3). Then, since (3.3) is used every time, we 
will not recall this frequently. 



3.2 The BTZ solution 

As an example, we calculate the mass of the BTZ BH [10] with the metric (2.9). We take the 
Einstein parts of each of the superpotentials (A.l), (A. 5) and (A. 7), and, keeping in mind a 
3D consideration, calculate their 01-components 
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(3.6) 



where the prime means d/dr. Taking into account a background with / = — r 2 A , for which 
f — f = —fi, and substituting (3.4) - (3.6) into (3.2), we obtain, as r — > oo, the unique result 

M = I E T^ l d<\> = ^ , (3.7) 
6 



which is quite acceptable for the global mass of the BTZ BH (see, e.g., [13]). The canonical 
superpotential has already been checked for calculating (3.7) in [14], for the other superpo- 
tentials the result (3.7) could be considered as a nice test. Using the AdS background with 
/ = — r 2 A + 1 one obtains M = n(fx + 1)/%. 

3.3 The static KK solution 

Now let us turn to (2.6); since it is the solution of the EGB theory one should try to 
calculate the mass with using the full formulae (A.l), (A. 5) and (A. 7) for this theory. The 
full background metric is to be chosen as = g AB x r-Q^ ab . Many formulae below take place 
for arbitrary / in (2.6), although in specific calculations we choose / = r 2 // 2 + q/r — q/r+. 
Let us turn to the (D — 2)-dimensional surface integral (3.2). Really, the distant surface is 
considered in {d = 3)-dimensional spacetime only, whereas the integral over the (D — d = 3)- 
dimensional compact space could be interpreted as a constant, which 'normalizes' the 6D 
Einstein constant Kq to the 3D one k 3 . Indeed, one has for the global mass constructed by 
(3.2): 

M=l dx D - 2 f^~ D T^ = £ d4>f^ d l% I dx D - d J-g D _ d = V ro <f tyypfiT™. 

JdT, v Jr^oo v Jro v Jr~^oo v 

(3.8) 

Thus, since T ^ ~ l//c 6 one could set k 3 = K§/V ro . At first we follow this prescription. 

With our assumptions, we find out that the Einstein parts of the 01-components of the 
superpotentials (A.l), (A. 5) and (A. 7) for the solution (2.6) are described only by the d- 
sector. Therefore, to calculate the Einstein parts, it is sufficient to use Eqs. (3.4) - (3.6), 
but only with y/—g 3 / k 3 replaced by yj— go/ K 6 ■ For all cases, in the natural background, 
the Einstein part in (3.8) gives a result corresponding to (3.7): 

M E = irjlV ro /K 6 . (3.9) 

We now construct the GB 01-components of the superpotentials (A.l), (A. 5) and (A. 7) 
for the solution (2.6). They consist of two parts. The first one is pure (d = 3)-dimensional: 

(3.10) 

m-f-n, (3.H) 

(3.12) 

(for brevity we suppress the subscript 'GB'). For / = r 2 /l 2 +q/r—q/r + , the behavior of (3.11) 
as r — > oo is ~ 1/r 3 , thus each of the variants (3.10) - (3.12) gives a zero contribution into 
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the integral (3.8). The other part of the GB 01-components is determined by the intersecting 
terms of the (d = 3)-sector and the scalar curvature of the (D — d = 3)-sector (2.8): 
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(3.13) 
(3.14) 
(3.15) 



where the subscript l (D~d) means that a quantity is without pure '-terms. We remark 
that both (3.10) and (3.13) are unique for each of (A. 3) and (A. 4). The asymptotic of each 
of (3.13) - (3.15) at spatial infinity in the natural background is ~ — jl, and their substitution 
into (3.8) gives the unique result: 



M GB = -TT[iV ro /KQ . 



(3.16) 



Thus, keeping in mind (3.9) one can see that the global mass defined in the natural back- 
ground by the total integral (3.8) is zero in all the three approaches. The same result is valid 
if the AdS background with / = r 2 // 2 + 1 is chosen. 1 

At least, this result could be anticipated for the field-theoretical approach. Indeed, the 
superpotential (A. 7) can be connected directly with the linearized equations [7]. Contracting 
the latter with £ a in (3.3), one selects the ci-sector only. However, under the present assump- 
tions, the tensor in (2.3) is equal to zero identically, therefore its linearization is equal to zero 
identically as well. This conclusion is supported by combining the expressions (3.6), with 
the replacement y/—g 3 /n3 — > y/^Sr)/ K 6, (3-12) and (3.15), which leads to zero identically. 
At the same time, the canonical and Belinfante corrected approaches give a zero result only 
asymptotically. 

Of course, the zero result cannot be acceptable. Analyzing (2.6), one can find out that, 
considering this system from the point of view of the Newtonian-like limit in 3 dimensions 
(see, e.g., [13]), this system must have a total mass. Thus one should conclude that a vacuum 
6D interpretation (2.2) with (3.8) is not successful. By this argument, one should consider the 
3D Einstein interpretation (2.12) with a created 'matter'. Calculating the global conserved 
quantity basing on (3.2), we can use only the surface integral, whereas a source (maybe not 
determined explicitly, as in (2.12)) is included into the current in the volume integral. Thus, 



1 The zero result has been recently obtained for a similar situation by other methods as well by R.G. Cai, 
L.M. Cao, and N. Ohta, "Black holes without mass and entropy in Lovelocj gravity", Phys. Rev. D, 81, 
024018; (Preprint arXiv:0911.0245 [hep-th]). 
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considering the solution (2.6), we can be restricted to only the Einstein parts of each of 
the superpotentials (A.l), (A. 5) and (A. 7) related to the non-vacuum equations (2.12). As 
a full background metric, one must again consider g AB in (2.6) (without r^ab); we choose 
/ = r 2 // 2 + q/r — q/r + again and use the Killing vector (3.3). Then, since the parameter q 
describes a 'created matter' in (2.12), such a background is not vacuum in 3 dimensions now. 
Nevertheless, the meaning of the notion 'real vacuum' is not changed, although it could be 
called wider as a 'real background' now. Also, the applied formalism remains powerful in 
non-vacuum backgrounds, and the structure of the superpotentials remains the same. Then 
again we use (3.4) - (3.6) and obtain the acceptable result of the type (3.7): 

M = njl/K 3 . (3.17) 

If AdS space with / = r 2 // 2 + 1 is chosen as a background, the mass of the system is 
M = 7r(/x + 1)/k 3 . Note that in both cases the parameter q makes no contribution. 

3.4 The radiative Vaidya KK solution 

For the radiative solution (2.10) we have carried out calculations similar to those in Sub- 
section 3.3. Though, in this case the lightlike f-coordinate is used instead of the time 
t-coordinate. We again calculate 01-components for the superpotentials, however, now E in 
(3.2) is defined as x° = v = constant, and the mass calculation is related to null infinity. In 
Eqs. (3.18) - (3.23) below, an arbitrary / = /(r) is considered. However, now there is no 
sense to connect a background (which must be static) with a horizon (which is changed in 
time). Therefore, in specific calculations we consider the AdS background with / = r 2 // 2 + l 
only. 

We first derive out the Einstein parts of all superpotentials: 

E I% = E 7% = E Xf = -^(/ - 7) (3.18) 

where / = f(v,r), which looks surprisingly simple, see, e.g., (3.4) - (3.6). The GB 01- 
components of the superpotentials (A.l), (A. 5) and (A. 7) for the solution (2.10) consist of 
two parts again. The pure (d = 3)-dimensional part is 
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For the AdS background one has as r — > oo: for (3.19) ~ 1/r 3 and for (3.21) ~ 1/V 2 , thus 
all (3.19) - (3.22) again give a zero contribution into the integral (3.8). As in the static 
case, the other part of the GB 01-components is determined by the intersection terms of the 
(d = 3)-sector and the scalar curvature of the (D — d = 3)-sector (2.8): 

-fOl -7-01 -7-01 V~9d ( f -J\ oq\ 

(D-d)-L C = {D-d)-L B = {D-d)-L s = ~2kqt~ ~ ' ' ^ ' 

One can see that these components precisely compensate the components (3.18). Thus, as 
in the previous subsection, the global mass defined in 6 dimensions is zero. Then one should 
follow the interpretation of the static case and reject the vacuum 6D derivation (2.2) with 
(3.8) as unacceptable one. 

We again consider Eq. (2.12) as a governing one. Restricting ourselves to the (i-sector only 
and repeating the steps of Subsection 3.3, we obtain in the AdS background M = 7r(/z(t> ) + 
l)//t3. This is in a correspondence with the static case. The mass flux for the radiating metric 
(2.10) is obtained simply by differentiating with respect to v : M = 7r/i(f)//t3. Comparing 
with the known BMS flux derivation [15], this looks acceptable. 

Concluding the section we assert that since the KK BH objects have classically defined 
global mass and flux, they bring 'matter' created by extra dimensions and a special structure 
of the objects themselves. If we set q = and q(v ) = 0, then, at least in the static case, 
Tab = in (2.2). However, this does not influence on our assertion because in all the cases 
q and q{v ) do not contribute into the global mass. Thus, mass/matter is created in a more 
wide sense than creating Tab in (2.2). 

4 Concluding remarks 

We will first discuss a well-known ambiguity in the canonical approach related to a choice 
of a divergence in the Lagrangian. We consider this problem in [9] and do not make a 
definite choice between [14] (or (A. 3)) and [16] (or (A. 4)). Indeed, both choices give an 
acceptable mass for the Schwarzschild-AdS BH tested in [9]. Here, the study of KK objects 
also does not give an answer because in all cases we have a unique result. However, in 
[16] arguments in favor (A. 4) are given. In multitimendional GR, the Katz and Livshits 
superpotential [16] turns out uniquely the KBL superpotential [17]; in EGB gravity, their 
superpotential naturally transfers into the KBL superpotential for D — 4. This is in a 
correspondence with the Olea arguments [18] where GB terms in the Lagrangian regularize 
conserved quantities even if D < 4. Lastly, the choice (A. 4) looks more preferable because 
(a) it is more 'symmetric' than (A. 3), (b) the canonical superpotential with (3.20) gives a 
zero global integral in 6 dimensions identically, as in the field-theoretical approach. 
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Now we turn to cosmological problems. As well known, the properties of dark energy and 
of dark matter are very weakly constrained by the cosmological observable data, therefore 
their derivation remains very uncertain. Thus a search for acceptable models describing the 
cosmic ingredients is very important, it is carried out very intensively, and even dramatically, 
see, e.g., the recent papers, reviews [19] - [27] and references there in. 

As an example, in the recent paper [28], recalling the 't Hooft ideas of 1985, so-called 
'quantum black holes' are discussed as elementary particles playing the role of the dark 
matter particles. The latter are assumed as weak interacting matter particles (WIMPs), 
which can have desirable TeV energies (see the aforementioned reviews). 'Quantum black 
holes' can be presented just like WIMPs, they can be stable and do not radiate in the 
Hawking-Bekenstein regime, unlike usual black holes. 

Our main results show that the solutions (2.6) and (2.10) have a classically defined mass 
and mass flux. This just presents a possibility for the KK BHs to be presented in the regime 
of 'quantum black holes'. Thus, the topic of the present paper, as we think, could be related 
to the dark matter problems. Concerning this, we remark the following. First, since the 
parameter q can describe additional (to gravity) interactions, its presence can suppress the 
WIMP idea. Then one needs to set q — 0, which is permissible, as has been remarked above. 
Second, it is desirable to have a positive mass for WIMP objects. We support this because, 
if a BH exists, one has fx > that leads to M > 0. Third, basing on the radiating regime, 
in [1] - [4] a scenario of forming KK BHs in EGB gravity was suggested. One could try to 
develop this scenario for various epochs. Keeping in mind all that, in future studies we plan 
an examination of more realistic models presented in [1] - [4]: they are 4D KK objects in 6 
and more dimensions of EGB gravity. 
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A Superpotentials in the EGB gravity 

In this Appendix, we represent an explicit form of the three types of superpotentials for 
perturbations in the EGB gravity [9]. The background quantities: Christoffel symbols r , 
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covariant derivatives D a , the Riemannian tensor R a Tp7r and its contractions are constructed 
on the basis of a background D-dimensional spacetime metric g pu . It is a known (fixed) 
solution of EGB gravity; the bar means that a quantity is a background one. One can find 
a detail derivation in [9] . We first present the superpotential in the canonical prescription: 



X a J = + GB X a J = K' 1 (g^D p e ] + g p[a A® r - Dl^]+e a E~ 

+ GB^ - GbW + K-^GBdP 

\g ap {prfpu + D v g pit - D p g uu ) and 2 



(A.l) 



where A a — T a — T 



~a/3 
GB% C 



K 

2aJ-g 



{r x ^ + 4g x[a R$ + 5 [ ?gfl x R\ D x ^ . 



(A.2) 



The vector density d x = e<1 x + GBd x could be defined as in [14] or following the prescription 
of [16]: 

d x = 2y^A%g x ^ + 4a^(R a ^ x -4R^g x ^ + 5^g x ^R)A^, (A.3) 
d x = 2y/=g~A [ ^g x] P + 4aJ=g~ (R a apx - 2R^g x ^ - 25 l «R x](3 + 5^g x]p R) . (A.4) 

The Einstein part in (A.l) is precisely the KBL superpotential [14, 17], which in 4D general 
relativity (GR) for the Minkowski background in the Cartesian coordinates and with the 



translation Killing vectors £ c 



is just the well-known Freud superpotential 



The Belinfante corrected superpotential in EGB gravity is 

xf = E xf + GB xf = kt 1 (e a Dj® x - D [a l^ + &°D x £n) + GB ri - gbW (A.5) 



where / a/3 = g a/3 — g a and 

GB r£ = -Dx {Ra Xa P + 4g x[a RP) + \2R T px[a - 2R px T [a - 8R x g p[a 



+ AR p g x ^ + Ag px R^ + 2R (S x g^ a - %g* a )] g® r g pa } £ 



2a 



K 



{R a XaP + 4g x[a R® + 8^g^ x k) D X C ■ 



(A.6) 



The Einstein part, e^e 13 , being constructed in arbitrary D dimensions, has precisely the form 
of the Belinfante corrected superpotential in 4D GR [30]. In the Minkowski background in 



2 Thc expression (A.2) is differed from the correspondent one in [9], where the mistake has been found. 
Nevertheless, the main results and conclusions in [9] are not changed; see Corrigendum: Class. Quantum 
Grav. 27 (2010) 069801 (2pp); Preprint arXiv:0905.3622 [gr-qc] . 
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the Cartesian coordinates and with the translation Killing vectors e&b ■> h transforms to the 
well-known Papapetrou superpotential [31]. 

Lastly, the superpotential in the field-theoretical derivation in EGB gravity is 

if = E lf + GB lf = k- 1 (^ v D [a h^ - $ a D v h® v + £ a D% - h v &D v Z® + \hD [a ^ 

+ | (2^D x NSS m ~ N a Js mX D x ^) . (A.7) 

where h a p = y/^(g a p - g aj3 ) and 

N p ^ p]v = (A.8) 

- {K \g u[X T ]p R + 2g p[x R p]u - 2g u[x R p]p - TT^] + (g p H^ - g u[x h p]p ) R 
+ 2 (h u[x W ]p - h p[x R p]u ) + 2 (h a[x g Ap R v a - h alx g^R p ) + 2 (h ap g^ x R$ - h al/ g p[x Wj) 

- 2g v[X g p]p K% + 4 (R a [M[p h v]c + R a [pu][x h p]a ) + 1h aT (iT^V 1 " - R c ^ x g^ v ) } . 

One obtains from (A.7) the Deser-Tekin superpotential [32] if one chooses the AdS back- 
ground. Again, doing simplifications in 4 dimensions as above, one obtains the Papapetrou 
superpotential [31] (note, see [8], that in 4D GR the Belinfante and field-theoretical ap- 
proaches give the same result). Under weaker restrictions, say, to AdS/dS backgrounds in 
4D GR, the superpotential (A.7) goes to the Abbott-Deser expression [33]. 
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